Stochastic control of quantum coherence 
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The application of a random modulation of a system parameter usually increases decoherence 
effects. Here we show how, employing an appropriate stochastic modulation, it is instead possible 
to preserve the quantum coherence of a system. 
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Controlling quantum coherence is one of the most fun- 
damental issues in modern information processing [Q. 
The most popular solution in the field of quantum infor- 
mation are quantum error correction codes ||] and error 
avoiding codes H , both based on encoding the state into 
carefully selected subspaces of a larger Hilbert space in- 
volving ancillary systems. The main limitation of these 
strategies for combatting decoherence is the large amount 
of extra space resources required Q ; in particular, if fault 
tolerant error correction is also considered, the number 
of ancillary qubits enormously increases. 

For this reason, other alternative approaches which do 
not require any ancillary resources have been pursued, 
and which may be divided into two main categories, ac- 
cording to the form of interaction with the system un- 
der study (5). If the interaction is one way, so that the 
controller acts on the system without obtaining any in- 
formation about its state, then the controller is called 
"open loop" ||. By contrast, if the controller acts on 
the system on the basis of information that it obtains 
about the state of the system, then it is called "closed 
loop" [^U^]. In standard open loop techniques, control of 
quantum dynamics is achieved through the application 
of suitably tailored, time-dependent and deterministic, 
driving forces. Here we want to extend open loop control 
strategies by considering the possibility of using stochas- 
tic parameter modulations for quantum control. 

The common wisdom is that whenever a system is 
subject to noise, the quality of the dynamic control is 
degraded, and that quantum coherence in particular is 
rapidly destroyed ||,[|. Here we show that this is not 
generally true and that quantum decoherence can be sig- 
nificantly suppressed if an appropriately tailored stochas- 
tic modulation of a system parameter is used. This fact is 
illustrated in this letter by considering the simple case of 
the dynamics of a single radiation mode in a lossy cavity, 
but the results can be generalized. In this open system, 
decoherence has a dissipative origin since it is due to pho- 
tons' leakage out of the cavity, and the stochastic control 
strategy will be implemented by modulating the cavity 
length. 

Let us consider a single radiation mode with annihila- 
tion operator a within a lossy cavity, whose characteristic 
frequency is oj — mrc/ 'L, with n an integer number, c the 



speed of light and L the cavity length. If photons' leak- 
age occurs through a partially transmitting mirror, the 
decay rate will be given by 7 = cT/2L, with T the mir- 
ror's transmittivity. 

In the case of optical frequencies, thermal excitation 
from the environment of the continuum of modes outside 
the cavity is negligible and the dynamics is well described 
by the master equation [|l 0f| 



p = Cp = —iuj [a 1 a, pi + ^T>[a]p . 



(1) 



where V[A]B = ABA'' - {A^A, B}/2 describes photon 
decay into the vacuum. This decay is also responsible 
for the rapid decay of any eventual quantum coherence 
generated within the cavity jnj . 

Let us now try to preserve the quantum coherence 
of the radiation mode using an appropriate stochastic 
control strategy. In particular, we randomly modulate 
the cavity length, that is, L — > L(t) = Lo/y(t) with 
y(t) a positive stochastic process. This is equivalent 
to a simultaneous random modulation of both the fre- 
quency and the decay rate of the cavity. This stochas- 
tic modulation of the cavity length moreover yields a 
dynamics which is indistinguishable from that driven 
by the constant, unmodulated, Liouvillian superoperator 
Co = —iujo [a' a, ■ ■ ■] + 'Jo'D [a] ■ ■ ■ , where the parameters 
ujq and 70 are fixed, in the presence of a random evolution 
time t' . In fact, one has 

p(t) = Texp|^ dftC(«)jp(0) =exp{£ o t'(t)}p(0), 

where T denotes time ordering, and we have defined the 
stochastic evolution time t'{t) — f Q dsy(s). This obser- 
vation reminds the recently proposed model-independent 
approach to decoherence in quantum mechanics in 
which the evolution time is regarded as a random vari- 
able. 

To establish a connection between the randomized time 
evolution of Refs. and the model of a cavity mode 

with a stochastically modulated cavity length we assume 
that the statistical properties of the cavity length mod- 
ulation factor y(t) are determined just by the Gamma 
probability distribution P(t, t') for the random evolution 
time t' of Ref. |^,|l2|- To be more specific, we assume 
that at discrete times t n separated by a time interval At, 
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independent random variables y(t n ) are generated, e.g. 
by a computer, according to the distribution 



P(y) 



¥) 



At/r yAt/r-1 



-yAt/r 



(2) 



T(At/r) 

The random number y(t n ) determines the "instanta- 
neous" cavity length L(t n ) = Lo/y(t„). Equation (||) 
is a Gamma probability distribution [fl3| , with parame- 
ter At/r, where t quantifies the strength of the fluctu- 
ations. In fact it is (y(t n )) = 1 and ((y(t n ) — l) 2 ) = 
r/At. Choosing the probability distribution (^|) means 
choosing a specific, uncommon way of modulating the 
cavity length. In fact it is easy to see pjjil^ that 
it implies a strongly non-Gaussian cavity length mod- 
ulation, which assumes Gaussian properties {P{y) — 
exp [-{y - l) 2 At/2r] /y/2wT/At) in the limit At/r > 1 
only. 

The above introduced discreteness in time, is dictated 
by the unavoidably finite rate of random number gener- 
ation. Nevertheless, the time interval At can be much 
smaller than the typical time scale upon which one ob- 
serves the system dynamics (this is essentially deter- 
mined by Jq 1 times the inverse of the mean photon num- 
ber). Hence, we can consider the independent cavity 
length rescalings y(t n ) as occurring continously in time, 
i.e. At — > 0. In the continuous approximation, y(t) be- 
comes a white, non- Gaussian stochastic process, which 
can be rewritten as y(t) = 1 + £(t), where is a zero- 
mean stochastic process, such that (£(i)£(t')) = TS(t—t'). 

As noted above, the evolution in the presence of the 
cavity length modulation can be reinterpreted as the dy- 
namics of a cavity with fixed length and with a random 
evolution time t'(t). However, since the sum (integral) 
of independent Gamma-distributed processes is again a 
Gamma-distributed process with the parameter given by 
the sum (integral) of the single parameters , the prob- 
ability distribution of the effective random time t'(t) is 
just the Gamma distribution of Refs. |p|,|l2|, 



p(t,t') = e - [tlT \ t>>0, 

r r(t/r) 



(3) 



with parameter t/r. This fact simplifies the study of the 
dynamics of the dissipative radiation mode in the pres- 
ence of the stochastic modulation of the cavity length. 
In fact any dynamical quantity can be obtained by first 
evaluating it in the absence of modulation and then av- 
eraging it over the random time distribution ([}]). 

A first interesting quantity is the time evolution of 
the intracavity mean photon number (a^(t)a(t)) = (n(t)} 
(the time average is denoted by the overbar). In the 
absence of any stochastic modulation one has (n(t)) = 
(n(t)} = (n(0)} exp(— 7 i), showing the energy decay due 
to photon leakage through the partially transmitting mir- 
ror. In the presence of the cavity length modulation one 
has instead 
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with the new effective energy decay rate j e given by 

le = r" 1 log(l + 7o r) . (4) 

It is always "f e < 70 and therefore cavity length modula- 
tion always yields inhibition of dissipation. What is rele- 
vant is that the suppression of energy damping increases 
for increasing fluctuation strength parameter r and that 
one has perfect inhibition, j e = 0, in the limit jqt — > 00, 
when the modulation becomes strongly non-Gaussian. 

Another interesting quantity is the behavior of the cav- 
ity fie ld, which is essentially expressed by the average 
(a(t)). In the absence of any stochastic modulation one 

has (a(t)) = (a(t)) = (<z(0)) exp(-iw * - 7o*/2). This 
quantity shows the effect of photon leakage on the phase 
of the intracavity field and it may provide some infor- 
mation on the possibility to control decoherence using 
cavity length modulation. In fact, since in this model 
decoherence is just caused by photon leakage, we expect 
that any control exerted on the field decay rate will re- 
flect itself into a control of quantum decoherence. In 
the p resence of the cavity length modulation one has 
(a{t)) = (a(0))e _<o * _;?t / 2 , where the renormalized field 
decay rate 7 is 



= r- 1 log (l + 7o r/2) 2 + 
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and the renormalized oscillation frequency is Co = 
t -1 arctan [wor/(l + 7or/2)]. The effective field decay 
rate 7 of Eq. (^) is different from the effective energy 
decay rate j e because it is sensitive not only to the de- 
cay rate modulation induced by the cavity length mod- 
ulation (as 7 e ), but also to the simultaneously induced 
frequency modulation, which may provide phase damp- 
ing effects. The behavior of the effective field decay rate 
(H) , is plotted in Fig. 1 as a function of the modulation 
strength parameter jqt. Curve (a) of Fig. 1 refers to 
Eq. (5) and, at variance with what happens for j e of 
Eq. (4), shows an initial increase of the effective cavity 
decay rate for increasing modulation amplitude r. This 
means that for not too large r, the modulation of the 
cavity length increases the cavity field decay rate. This 
decay acceleration reaches a maximum around ojqt ~ 1 
and then starts to decrease for increasing r. What is 
interesting is that the ratio 7/70 becomes less than one 
and tends to zero for larger t. This means that not only 
dissipation, but also field decay can be completely in- 
hibited by the cavity length modulation at a sufficiently 
large r parameter, when the stochastic modulation as- 
sumes strongly non-Gaussian properties. The threshold 
value Tth for decay inhibition, 7 < 70, depends on the 
cavity quality factor Q = luq/^/q, and it can be calculated 
iteratively, obtaining 



dt' P(t,t')(n(t')) = (n(0))e 



-7e* 



7o r t/l «21ogQ + 21og [21ogQ] 



(6) 
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Curve (b) of Fig. 1 instead refers to the cavity field decay 
rate in the case of a Gaussian cavity length modulation, 
whose expression essentially coincide with that obtained 
by extrapolating the expansion of Eq. (^|) at first order 



in r, that is, 7c a 



7o 



7o/4)r. The Gauss lan 



modulation case always yields an accelerated decay rate 
(it is always ujq 3> 70 in optical cavities) showing that 
the statistical properties of the stochastic modulation (its 
non-Gaussian properties in particular) are a fundamental 
ingredient for achieving a significant decay inhibition. 
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FIG. 1. Log-log plot of the ratio 7/70 as function of 70T. 
Curve (a) refers to the Gamma stochastic modulation of 
Eq. (|^), and curve (b) refers to a Gaussian stochastic cav- 
ity length modulation. We have also set Q = 10 2 . 

Let us now directly address the decoherence control 
issue. We consider as initial state of the cavity field a 
linear superposition state. In order to control the coher- 
ence in the continuous variable case we shall consider the 
well known Schrodinger cat state, a superposition of two 
coherent states of the form \a) + 1 — a) and see what hap- 
pens by employing the above stochastically modulated 
dynamics. The same could eventually be done for a su- 
perposition of Fock states. 

The time evolution of the Schrodinger cat state in 
the absence of any modulation is determined by the 
usual Liouvillian and it can be described in the following 
way m P(t) = N- 2 {[\a(t)){a(t)\ + | - a(t)){-a(t)\] + 
exp[-2T«| 2 (l ~v(t))][Ht))(-a(t)\ + \ - a(t))(a(t)\]} , 
where we have introduced a(t) = aexp[— (iujq + 7 /2)i] 



and 77 (i) = e 



-70* 



A good characterization of the time 



development of the quantum coherence of the state of 
the cavity mode is provided by the visibility with re- 
spect to an observable For the quadrature observ- 
able X = (a + a))/^/2 1 the quantum visibility is given by 



v 



e -2| Q | 2 [i-^(t)] (x\a(t)){-a{t)\X) 



^[(X\a(t))(a(t)\X)} [{X\ a(t))(-a(t)\X)} 



(7) 



where 
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In the absence of modulation, Eq. (]?]) leads to the sim- 
ple result V = exp{— 2\a\ 2 [1 — I n the pres- 
ence of the stochastic modulation of the cavity length, 
the corresponding visibility can be evaluated by per- 
forming an appropriate average of the dynamical quan- 
tities over the probability distribution P(t,t'). In 
particular, we have to consider the following replace- 
ments in Eq. (fl) e -2\ a \ a ^-vW ) (x\g{t)) (-ct(t)\X) 
e-2\<*\ 2 (i-v(t))(X\ a(t))(-g(t)\X), (X\a(t))(a{t)\X) -» 
(X\a(t))(a(t)\X), (X\ - a(t)){-a(t)\X) 
(X\ — a(t))(— a(t)\X) to get the corresponding, aver- 
aged, V. A cumbersome analytic expression can be 
obtained jTJJ] and the corresponding behavior of V as a 
function of time for different values of the modulation 
strength parameter t is shown in Fig. 2. The relevant 
result is that the visibility, i.e., the quantum coherence 
properties of the system, behaves in the same way as the 
field decay rate. In particular we see either an accelera- 
tion, or, more importantly, a deceleration of decoherence 
according to the value of the parameter r. The usual de- 
cay of the visibility in the absence of modulation (r = 0) 
is shown with a dashed curve. When j t ^ we ob- 
serve an acceleration of the decay of the visibility (lower 
curve) when the modulation strength r is not too large 
(70T = 1.5 in the figure) or a slowing down of the de- 
cay (upper curves) when r becomes sufficiently large 
(7o T = 20, 100 in Fig. 2). The threshold value between 
the two behaviors coincides with that for decay inhibition 
T th ofEq. (§). 




FIG. 2. Time evolution of the visibility for different val- 
ues of 70T. Curves from bottom to top refer to 70T = 1.5, 



7or = 0, 70T = 20, 70T 
a = 2i. 



10 respectively. We have also used 



The above results show that cavity dissipation, field 
decay, and decoherence can be inhibited if an appro- 
priate random modulation of the cavity length is ap- 
plied. We have assumed a Gamma distributed stochastic 
process (see Eq. (g)) and we have seen that the inhi- 
bition becomes significant when the stochastic modula- 
tion becomes strongly non-Gaussian (70T large). Since 
the statistical nature of the modulation plays such a rel- 
evant role, it is important to establish the stability of 
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the above decoherence control results with respect to 
small changes of the stochastic modulation. As suggested 
above, the Gamma stochastic modulation of the cavity 
length could be experimentally implemented using a cav- 
ity with a computer-controlled length, and a fast ran- 
dom number generator. Any eventual imperfect stochas- 
tic modulation of the cavity length can be described in 
terms of an additional, zero-mean, white, stochastic pro- 
cess e(i), such that (e(i)e(t')) = aS(t — t'). The stochas- 
tic process e(t) describes the "error" in the modulation 
at time t, so that the actual fluctuating cavity length 
is L(t) = Lo/(y(t) + e(t)). We have checked that the 
above decay and decoherence inhibition results are not 
significantly changed as long as the strength u of the 
additional, undesired, modulation e(t) is not too large. 
More precisely, it is possible to see with a perturbative 
treatment that the results are stable if a^foQ 2 <C 1. This 
condition can be explicitely seen for example in the case 
of the renormalized field decay rate 7 and it is shown in 
Fig. 3 in the case of a Gaussian stochastic process e(t). 
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FIG. 3. Log-log plot of the ratio 7/70 as function of 70T 
for different values of the "error strength" a. Curves are for 
(770 = (a), 10" 5 (b), 10" 4 (c). We have also set Q = 10 2 . 

In conclusion, we have studied the possibility of a 
stochastic control of (dissipative) decoherence by tailor- 
ing suitable random modulations of a system parameter. 
Against the widespread opinion that "noise" is detrimen- 
tal for quantum effects, we have shown that if the sta- 
tistical properties of the modulation are appropriately 
chosen, this stochastic control strategy could be used in 
principle to control decoherence. Here we have consid- 
ered the specific model of a single cavity mode with a 
randomly modulated cavity length. We have seen that, 
when the modulation is stochastic, with strongly non- 
Gaussian properties, decoherence and dissipation can be 
inhibited and that the scheme can even tolerate some 
imperfection in the modulation. Although we have con- 
sidered a specific model, our results can be generalized 
to dissipative system in Ohmic environments where the 
damping rate and the frequency have the same fluctua- 
tions |l0|. That allows us to recast the above described 



treatment. Finally, our approach shares some similarities 
with the inhibition of atomic decay through random ac- 
Stark shift discussed in Ref. However our proposal 
is different since it strongly depends on the statistical 
properties of the random modulation and it is especially 
suited to the control of quantum decoherence. Another 
analogy occurs with the use of kicks to prevent the decay 
of a system Q. In this latter case, dephasing introduced 
by kicks were deterministic processes well defined in time. 
Instead, the present approach is merely probabilistic, so 
it would be more manageable. 
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